Flavor Changing Neutral Currents Transition of the Eg to Nucleon in Full QCD and 

Heavy Quark Effective Theory 



K. Azizii't, M. Bayar^'t, M. T. Zeyrek^-*, 
Physics Division, Faculty of Arts and Sciences, Dogu§ University, Acibadem-Kadikoy, 

34722 Istanbul, Turkey 
^ Department of Physics, Kocaeli University, 41380 Izmit, Turkey 
* Physics Department, Middle East Technical University, 06531, Ankara, Turkey 

^kazizi@dogus. edu.tr 
^melahat. bayar@kocaeli.edu.tr 
•^zeyrekQmetu. edu.tr 

The loop level flavor changing neutral currents transitions of the Ej, — > n and T,c ^ p 
are investigated in full QCD and heavy quark effective theory in the light cone QCD sum rules 
approach. Using the most general form of the interpolating current for Eg , Q = 6 or c, as members 
of the recently discovered sextet heavy baryons with spin 1/2 and containing one heavy quark, the 
transition form factors are calculated using two sets of input parameters entering the nucleon distri- 
bution amplitudes. The obtained results are used to estimate the decay rates of the corresponding 
transitions. Since such type transitions occurred at loop level in the standard model, they can be 
considered as good candidates to search for the new physics effects beyond the SM. 
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I. INTRODUCTION 



The — n l^l^ and 'Sc P are governed by flavor changing neutral currents (FCNC) transitions of 6 — d 
and c — ?► w, respectively. These transitions are described via electroweak penguin and weak box diagrams in the 
standard model (SM) and they are sensitive to new physics contributing to penguin operators. Looking for SUSY 
particles [l|, light dark matter [2] and also probable fourth generation of the quarks is possible by investigating such 
loop level transitions. This transitions are also good framework to reliable determination of the Vtb, Vtd, Vcb, and Vhu 
as members of the Cabibbo-Kobayashi-Maskawa (CKM) matrix, CP and T violations and polarization asymmetries. 
The Sf, c as members of the spin 1/2 sextet heavy baryons containing a single heavy bottom or charm quark are 
considered by their most general interpolating currents which generalize the loffe current for these baryons. In the 
recent years, important experimental progresses has been made in the spectroscopy of the heavy baryons containing 
heavy b or c quark P-[l0|. Having the heavy quark makes these states be experimentally narrow, so their isolation 
and detection are easy comparing with the light baryons. Experimentally, investigation of the semileptonic decays of 
the heavy baryons, may be considered at large hadron collider (LHC) in the future, hence theoretical calculations of 
the decay properties can play crucial role in this respect. 

In our two recent works, we analyzed the tree level semileptonic decays of Sb to proton pH and Ah(Ac) — ?> p{n)lv 
in light cone QCD sum rules. In full theory, these tree level transitions in the SM are described via six form 
factors (for details and more about the works devoted to the semileptonic decays of the heavy baryons using different 
phenomenological methods see [ll|, [l2| and references therein). In the present work, considering the long and short 
distance effects, we calculate the 12 form factors entering the semileptonic loop level — n l^l^ and Ec — > P ^^^^ 
transitions using the light cone QCD sum rules in full theory as well as heavy quark effective theory (HQET). The 
short distance effects are calculated using the perturbation theory and long distance contributions are expanded in 
terms of the nucleon distribution amplitudes (DA's) with increasing twists near the light cone, ~ 0. We use 
the value of the eight independent parameters entering to the nucleon DA's from two different sources: predicted 
using a simple model in which the deviation from the asymptotic DAs is taken to be 1/3 of that suggested by the 
QCD sum rule estimates [l^ and obtained via lattice QCD |l3 - [T6j . Using the obtained form factors, we predict the 
corresponding transition rates. Investigation of these decays can also give essential information about the internal 
structure of Yj^^c baryons as well as the nucleon DA's. 

The layout of the paper is as follows: in section II, we introduce the theoretical framework to calculate the form 
factors in light cone QCD sum rules method in full theory. The HQET relations among the form factors are also 
introduced in this section. Section III is devoted to the numerical analysis of the form factors and their extrapolation 
in terms of the transferred momentum squired, , their HQET limit and our predictions for the decay rates obtained 
in two different sets of parameters entering the nucleon distribution amplitudes. 
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II. LIGHT CONE QCD SUM RULES FOR TRANSITION FORM FACTORS 

Form factors play essential role in analyzing the ^ n and T,c ^ p transitions. At quark level, these 
decays proceed by loop 6 (c) — )• (u) transition, which can be described by the following effective Hamiltonian: 



^eff 



Gf aVq'Q Vq^q 

2\/2 TT 



2mQ qia^^q^^l + J5)QIY 



(1) 



where, Q' refers to the u, c, t for bottom case and d, s ,b for charm case, respectively. The main contributions come 
from the heavy quarks, so we will consider Q' = t and Q' = b respectively for the S(, ^ n Z+Z~ and Sc ^ p tran- 
sitions. The amplitude of the considered transitions can be obtained by sandwiching the above Hamiltonian between 
the initial and final states. To proceed, we need to know the the matrix elements (7V| J^'''-^|Sq) and {N\Jj^'^^\T,Q), 
where Jj['^{x) = q{x)'yij,{l — ^5)Q{x) and J^'^^{x) = q{x)iani,q"{l + 75)Q(a;) are transition currents entering to the 
Hamiltonian. From the general philosophy of the QCD sum rules, to obtain sum rules for the physical quantities we 
start considering the following correlation functions: 



n'^{p,q) = i j d^xe^'^Nip) I T{Jl^^\x)J^'^m \ 0), 
= i j d^^e'^^Nip) I T{4^''\x)J^'^m \ 0), 



(2) 



where, J^'^ is interpolating currents of Sb(c) baryon and p denotes the proton (neutron) momentum and q = (p + q) —p 
is the transferred momentum. The main idea in QCD sum rules is to calculate the aforementioned correlation functions 
in two different ways: 

• In theoretical side, the time ordering product of the initial state and transition current is expanded in terms 
of nucleon distribution amplitudes having different twists via the operator product expansion (OPE) at deep 
Euclidean region. By OPE the short and large distance effects are separated. The short distance contribution 
is calculated using the perturbation theory, while the long distance phenomena are parameterized in terms of 
nucleon DA's. 

• From phenomenological or physical side, they are calculated in terms of the hadronic parameters via saturating 
them with a tower of hadrons with the same quantum numbers as the interpolating currents. 

To get the sum rules for the physical quantities, the two above representations of the correlation functions are 
equated through the dispersion relation. To suppress the contribution of the higher states and continuum and isolate 
the ground state, the Borel transformation as well as continuum subtraction through quark-hadron duality assumption 
are applied to both sides of the sum rules expressions. 

The first task is to calculate the aforementioned correlation function from QCD side in deep Euclidean region where 
(p + q)'^ <C 0. To proceed, the explicit expression of the interpolating field of the Sq baryon is needed. Considering the 
quantum numbers, the most general form of interpolating current creating the Sg from the vacuum can be written 
as 



J^Ha;) = \^qI''{x)CQ''{x)^j,q'2{x)-^Q'^'^{x)CqUx)^j5qt{x) 

+13 ^{q^''{x)C^5Q''{x)}ql{x) - {Q^''{x)C^5qUx)}q1{x)^ , (3) 

where, C is the charge conjugation operator and 8 is an arbitrary parameter with 8 = —1 corresponding to the loffe 
current, qi and q2 are the u and d quarks, respectively and a, b, c are the color indices. Using the transition currents, 
and J^^and contracting out all quark pairs via the Wick's theorem, we obtain the following representations of the 
correlation functions in QCD side: 



-j=e'''"' I d'^xe"'''< 



7/^(1 - 75) 



' Sq {-x)0a {Nip)K {0)4 {xH (0) 1 0) , 



(4) 
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(C75 (^)p»? 



(C);3r,(75)p0 - (C)0/3(75)p,) + 



(C75)/3r,(/)p0 



ia^,q"{\ + 75) ^ 5Q(~x)^.(iV(p)M:;(0)4(a;)^i;(0)|0) 



where, Sq{x) is the heavy quark propagator and its expression is given as [l7| : 



(5) 



Sq{x) = S[j''"'{x) - ig. 



dv 



{ml - fc2)' 



U2 "-"M 



where 



(6) 



(7) 



and -fCi are the Bessel functions. When doing calculations, we neglect the terms proportional to the gluon field 
strength tensor because they are contributed mainly to the four and five particle distribution functions and expected 
to be very small in our case [Tsl-Eoll. The matrix element {N{p) \ e°*'^(i°(O)d^(x)u0(O) | 0) appearing in Eqs. (|4|5p 
denotes the nucleon wave function, which is given in terms of some calligraphic functions ^ISi, .18, - . 21.] : 



4(0|e'^^^<(aia:)4(a2x)u^(a3a:)|A^(p)) 

2 2 

+ VimN{-fr^C)^pN^ + V2m%{-fr^C)^p{^N)^ + (Vi + '^^^V^'){i>C)^p{-i5N)^ 

+ V2mN{i>C)ap{^-1^N)^ + V3TOAr(7MC)a/j(7^75A^)7 + V4™?^(^C)a/3 (75^^)7 

+ V5TOir(7MC)a/3(*^'"' 2:^75^)7 + V6m^(;^C)„;3 (^75^^)7 + Ml 

2 2 

+ ^^^f )(;^75C)a/jiV^ +-42mw(|^75C)a^(^iV)^ +^3m^v(7M75C)„^(7'^7V)^ 



(Ti 



+ r7m2^(a^^C)c«;3K'';^75^)7 +T8m^(xV^^C)„^(cr^%75A^)-y 



(8) 



The calligraphic functions have not definite twists but they can be expressed in terms of the nucleon distribution 
amplitudes (DA's) with definite and increasing twists by the help of the scalar product px and the parameters a^, 
i = 1,2,3. The explicit expressions for scalar, pseudo-scalar, vector, axial vector and tensor DA's for nucleons are 
given in Tables HI HIl Hill HVl and IVl respectively 



Si 



2pxS2 = Si - 5*2 



TABLE I: Relations between the calligraphic functions and nucleon scalar DA's. 
Each distribution amplitude F{aipx)= Si, Pi, Vi, Ai, Ti can be expressed as: 

F{aipx) = / dxidx2dxz5{xi + 2:2 + xa - l)e*P''^'^''''F(a;^) . 



(9) 
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Vi=Pi 

2pxV2 =Pl-P2 



TABLE II: Relations between the calligraphic functions and nucleon pseudo-scalar DA's. 



Vi=Vi 

2pxV2 = Vl-V2-V3 

2V3 = V3 

4paV4 = -2Vi + ^3 + V4 + 2Vt, 

4pa:V5 ^Va-Vs 

4(pa)^V6 = -1/1 + V2 + Va + VA + Vr,-V(, 



TABLE III: Relations between the calligraphic functions and nucleon vector DA's. 



where, Xi with z = 1, 2 and 3 are longitudinal momentum fractions carried by the participating quarks. Using the 
nucleon wave functions, which their explicit expressions are calculated in [isj and the expression for the heavy quark 
propagator, and after performing the Fourier transformation, the final expressions of the correlation functions for both 
vertexes are found in terms of the nucleon DA's in QCD or theoretical side. For simplicity, we present the explicit 
expressions of the nucleon DA's in the Appendix. 

The next step is to calculate the phenomenological or physical sides of the correlation functions. Saturating the 
correlation functions with a complete set of the initial state, isolating the ground state and performing the integral 
over X, we get: 

^i, , ^ {Nip) I J*--^(0) I EQ(p + g,,s))(EQ(p + g,5) I J^Q(O) | 0) 
TT/// ^ ^{N{p)\J';''\Q)\^Q{p + q,s)){^Q(p + q,s)\J^'^{0)\0) 

^ , -{p + ir 

where, the ... denotes the contribution of the higher states and continumii. The baryonic to the vacuum matrix 
element of the interpolating current, i.e., (Eq(p + q, s) | /^"(O) | 0) can be parameterized in terms of the residue, 
as: 

(EQ(p + g,s) I J^«(0) I Q) = \^^u^^{p + q,s). (12) 

To proceed, we also need to know the transition matrix elements, {N{p) \ Jjj^'^ \ T,q{p + q, s)) and {N{p) \ Jf^'^^ \ 
Sq(p + q,s)). In full theory, they are parameterized in terms of 12 transition form factors, fi, gi, and gj with 
i = 1 — > 3 by the following way: 

{N{p) I J';'\x) I Eq(p + q)) = N{p) [7^/1 (Q2) + za^,qV2(Q') + q'' h{Q^) - l^l,gi{Q^) ~ *^M-75g^92(Q') 

- q^l59s{Q'^) U^^{p + q), 

(13) 

and 

{N{p) I Jl:'"{x) I Sq(p + q)) = N{p) [7p/r(g') + ^<7^uq''f^{Q') + q^fliQ') + lf.l59l{Q') + ^'J^..l5q'' 9^ [Q') 

+ q'"759liQ^) u^Q{p + q), 

(14) 

where — ~q^. Here, N{p) and uy.q{p + q) are the spinors of nucleon and Eg, respectively. Using Eqs. pH)) . pT|) . 
([T2|). (|T3)) and (fT4)) and performing summation over spins of the Sq baryon using 

X! {p + g, ip + q, s) =i)+ i + , (15) 
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Ai = Ai 



2pxA2 



'Ai +A2~A3 



2A3 = A3 



4pxA4, = -2Ai -A3-A4 + 2A5 
ApxAs = A3~ A^ 



^{pxyAt, = Ai - yl2 + A3 + ^4 - ^5 + At 



TABLE IV: Relations between the calligraphic functions and nucleon axial vector DA's. 



2pxT2 = Ti + Ta - 2T3 



2T3 = T7 



2pxT4, = Ti - T2 - 2T7 



2px% = -Ti + T5 + 2rg 



MpxYTe = 2T2 - 2T3 - 2T4 + 2n + 2Tt + 2Ta 



ApxTr = Tt — Is 



AjpxyTs = -Ti +T2+T5 -n + 2T7 + 2Ts 



TABLE V: Relations between the calligraphic functions and nucleon tensor DA's. 



we obtain the following expressions 



^In - iP + 1)' 



(16) 



and 



As 



ra 



g"75ff3^(Q' 



{v + q)- 



:N{p) [iJ^iQ^) + i<j,,q^fUQ^) + q'^fliQ^) + 7^755^(0") + ^<J^.i5q''g^{Q^) 



4+ ^ + msg) + 



Using the relation 

Na^^q''uT,Q =i N[{mN + mT,Qhf_L- i2p + q)^]u^Q, 
in Eqs. (fT6| and Eqs. (fT7| . we attain the final expressions for the physical side of the correlation functions: 



(17) 



(18) 



ni(p,g) 



As. 



ip + q)' 



:N{p) 



2h{Q')p^ + - h{Q^){mN - msj + /2(Q')(m^ - m|j 7^ 



+ |/i(Q') - /2(Q')(mw + msQ)| 7p ^ + 2/2(Q2)p^ ^+ l^h{Q^) + h{Q^)^ (m^ + m^^)q^ 
/2(Q') + /3(0')|9m ^+25i(Q')Pm75- |5i(0')(mw+ms«) - 52(Q')(m^ - m|^)| 7^75 + 
9i{Q'^) - 92{Q'^){mM - tti^q) \"ff, ^75 + 252(Q^)Pm ^75 + \ 92{Q'^) + gsiQ'^) > {mN - ^2^)9^ 75 



+ 52(Q') + g3(Q') ^75 



(19) 
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and 



"^Eq - + 9) 



92m+9HQ )><l^^ ^75 



(20) 



In order to calculate the form factors or their combinations, /i, /2, f2 + fs, 9i, 92 and 92 +53, we will choose the 
independent structures p^, P;^^, q^i, Pm75' -Pm^75 5 a'^d q^^j5 from Eq. (jl9p . respectively. The same structures are 
chosen to calculate the form factors or their combinations labeled by T in the second correlation function in Eq. (|20p . 

Having computed both sides of the correlation functions, it is time to obtain the sum rules for the related form 
factors. Equating the coefficients of the corresponding structures from both sides of the correlation functions through 
the dispersion relations and applying Borel transformation with respect to (p + g)^ to suppress the contribution of 
the higher states and continuum, one can obtain sum rules for the form factors /i, /2, /s, 91, 92, 93, f\, /J, fj, 
91 1 92 a-iid 93- III heavy quark effective theory (HQET), where tuq — >■ 00, the number of independent form factors 
is reduced to two, namely, Fi and i^2- In this limit, the transition matrix element can be parameterized in terms of 
these two form factors in the following way [l^, [2^ : 

{N{p) I dVb I ^Q{p + q)) = iV(p)[Fi(Q2)+ jiF2iQ^)]Tu^,,{p + q), 

(21) 

where, F is any Dirac matrices and /i — -2^^. Here we should mention that the above relation is exact for A-like 



baryons, where the light degrees of freedom are spinless. For the E like baryons this relation cannot hold exactly and 
has to be replaced by a more complicated relation. In the present work, we will use the above approximate relation 
for the considered transitions. Comparing this matrix element and our definitions of the form factors in Eqs. (jl3p 
and (|14p. we get the following relations among the form factors in HQET limit (see also [2^ [25l|) 

/i = 51 = /2 = 52 = ^1 H F2 



F 



2 



/2 = 92 = /s = 93 
rT T ^2 2 

A = 5i = q 

13 = ('"Sq - mN) 

93 = — —{mT.Q+mN) (22) 

Looking at the above relations, we see that it is possible to write all form factors in terms of /i and /2, so we will 
present the explicit expressions for these two form factors in the Appendix and give extrapolation of the other form 
factors in finite mass as well as HQET in terms of in the numerical analysis section. 

The expressions of the sum rules for form factors show that we need to know also the residue As^ . This residue is 
determined in [26j: 

- Al^e^^^e^*'^^ = / e"Jp(s)ds + e^F, (23) 
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where, 



and 



{<dd> + <uu >) g4^2 j i;;^^^'^™ " -^^'^o^ - 6^-11 ) + 3mQ(2V'io - V'n " V'12 + 2i/'2i) 



20487r4 



[5 + /?(2 + 5/?)] [12^10 - 6V'20 + 2V'30 - 4V'4i + ^42 - 12Zn(^)], 



(24) 



24 



< ddxuu > 



2M| 4Af| 



(25) 



Here, ipn 



:™-(m^)" — 



are some dimensionless functions. 



III. NUMERICAL RESULTS 



This section deals with the numerical analysis of the form factors as well as the total decay rate of the loop level 
Sf, — > n£^£~ and — > p£~^£~ transitions in both full theory and HQET limit. In obtaining numerical values, 
we use the following inputs for masses and quark condensates: {uu){l GeV) — {dd){l GeV) — —(0.243)"^ GeV^, 
rrin = 0.939 GeV, nip = 0.938 GeV, rub = 4.7 GeV, = 1.23 GeV, m^, = 5.805 GeV, ms, = 2.4529 GeV and 
TO§(1 GeV) = (0.8 ± 0.2) GeV^ From the sum rules expressions for the form factors, it is clear that the nucleon 
DA's (see Appendix) are the main input parameters. These DA's contain eight independent parameters, namely, 
/at, Ai, A2, Vi, A^, ff, /" and /j- AH of these parameters have been calculated in the framework of the light 
cone QCD sum rules and most of them are now available in lattice QCD _14 ,16] (see Table rvT]) . Here, we should 
stress that in [l3| those parameters are obtained both as QCD sum rules and assymptotic sets, but to improve the 
agreement with experimental data on nucleon form factors, a set of parameters is obtained using a simple model in 
which the deviation from the asymptotic DAs is taken to be 1/3 of that suggested by the QCD sum rule estimates 
(see pL3]). We will use this set of parameters in this paper and refer it as setl (see Table rvT]) . In the following, we 
also will denote the lattice QCD input parameters by set2. 





setl 


set2 or Lattice QCD [14-16] 


/iV 


(5.0 ±0.5) X 10"=^ GeV^ 


(3.234 ±0.063 ±0.086) x 10"^ GeV^ 


Ai 


-(2.7 ±0.9) X 10-^ GeV^ 


(-3.557 ±0.065 ±0.136) x 10"^ GeV^ 


A2 


(5.4 ± 1.9) X 10"^ GeV^ 


(7.002 ±0.128 ±0.268) x 10"^ GeV^ 




0.30 


0.3015 ± 0.0032 ±0.0106 


AI 


0.13 


0.1013 ± 0.0081 ±0.0298 


ft 


0.33 




fi 


0.09 




h 


0.25 





TABLE VI: The values of the 8 independent parameters entering the nucleon DA's. The first errors in lattice values are 
statistical and the second errors correspond to the uncertainty due to the Chiral extrapolation and renormalization. For last 
tree parameters, the values are not available in lattice and we will use the setl values for both sets of data. 

The explicit expressions for the form factors also show their dependency to three auxiliary mathematical objects, 
namely, continuum threshold Sq, Borel mass parameter Mg and general parameter [3 entering to the most general 
form of the interpolating current of the initial state. The form factors as physical quantities should be independent 
of these parameters, hence we need to look for working regions for them. The working region for Borel mass squared 
is determined as follows: the upper limit of Mg is chosen demanding that the series of the light cone expansion 
with increasing twist should be convergent. The lower limit is determined from condition that the higher states and 
continuum contributions constitute a small fraction of total dispersion integral. Both conditions are satisfied in the 
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regions 15 GeV"^ < M% < 30 GeV'^ and 4 GeV'^ < M| < 12 GeV^ for bottom and charm cases, respectively. The 
value of the continuum threshold Sq is not completely arbitrary and it is correlated to the first exited state with 
quantum numbers of the initial particle interpolating current. Our numerical calculations show that the form factors 
weakly depend on the continuum threshold in the interval, (tosq +0.5)^ ^ sq < {rn^Q +0.7)^. To obtain the working 
region for /? at which the form factors are practically independent of it, we look for the variation of the form factors 
with respect to cos6 in the interval —1 < cosO < 1 which is equivalent to the — oo < /3 < oo, where /3 — tanO. As 
a result, the interval —0.5 < cos9 < 0.6 is obtained for /3 for both charm and bottom cases. In this interval, the 
dependency on this parameter is weak. 

The next step is to discuss the behaviour of the form factors in terms of the . The sum rules predictions for the 
form factors are not reliable in the whole physical region. To be able to extend the results for the form factors to 
the whole physical region, we look for a parametrization of the form factors such that in the reliable region which 
is approximately 1 GeV below the perturbative cut, the original form factors and their fit parametrization coincide 
each other. Our numerical results lead to the following extrapolation for the form factors in terms of q^: 



where the fit parameters a, b and mfu in full theory and HQET limit are given in Tables lVlIl IVIIIl IIXI and 1X1 using two 
sets for the independent parameters. These Tables, show poles of the form factors outside the allowed physical region. 
Therefore, the form factors are analytic in the full physical interval. In principle, we can use fit parametrization either 
with single pole or double poles. However, when we combine them the accuracy of the fitting becomes very high, 
specially when the pole is the same for two parts. We could start from fdq^) = 21 — j + t — 57 — j— 1 however for all 

form factors mfn gets too close to mi, so the fit becomes numerically unstable. In such a case, it is appropriate to 
expand the above relation to first order in mfn — mi, which gives the Eq. (I26p used to extrapolate the factors over 
the whole range of q^. For the same situation in B — > D mesonic transition see for instance [2^, [2^. The values 
of form factors aX q^ = are presented in Tables IXII and IXIII in both full theory and HQET for bottom and charm 
cases, respectively. In extraction of the values of form factors at q^ = 0, the mean values of the form factors obtained 
from the quoted ranges for the auxiliary parameters have been considered. When we look at these Tables, we see that 
although the values for the eight independent input parameters for two sets are close to each other but the results 
for the central values of some form factors differ in two sets, considerably. The numerical results show that the result 
of sum rules are very sensitive to these parameters specially /at, A2 and A". Within the errors, the quoted values 
become close to each other for both sets. The numerical analysis depicts also that all form factors approximately 
satisfy the HQET limit relations in Eq. ([22]) within the errors for both sets of input parameters and bottom case, 
Q = 6 at = 0. However for the charm case, Q — c although some of the relations are satisfied but most of them 
arc violated at = 0. This is an expected result since the nic — )■ 00 limit is not as reasonable as the nib 00. 

Our next task is to calculate the total decay rate of the FCNC Sfc — > p£^£^ and Sc — > nt^l" transitions in the 
full allowed physical region, namely, 4m^ < g^ < ("^s^ ^ ~ ™p.n)^- To derive the expression for the decay rate, we 
will make the following assumptions (see also [13]): the CLEO predicts the value R = ^ = —0.25 ± 0.14 ± 0.08 for 
the ratio of the form factors of Ac KeCi, at HQET limit [3l[. This result shows that IF2I < |i^i| and considering 
Eq. ((22)) . the form factors /i, gi, /i, /J and are expected to be large comparing to the other form factors since 
they are proportional to the Fi. Moreover, it is clear from the considered Hamiltonian as well as the definition of the 
transition matrix elements in terms of the form factors that the form factors labeled by T are related to the Wilson 
coefficient C7 which is about one order of magnitude smaller than the other coefficients entered to the Hamiltonian, 
i.e., Cg and Cio, hence their effects expected to be small. As a result of the above procedure, the following results for 
decay width describing such transitions is obtained [30| : 



— (Eg ^ Nl+l ) = ^^^-3 m^^ ^<i> {s)J 1 - —f (s) , (27) 



where 



Rsq (s) = ri (s) + r2 (s) + r3 (s) 



(28) 





setl 


set2 




a 


b 


rrifit 


a 


b 


rrifit 


fl 


-0.16 


0.29 


5.70 


0.027 


0.044 


5.02 




0.008 


-0.02 


5.96 


0.018 


-0.024 


7.96 


h 


0.011 


-0.024 


6.34 


-0.003 


-0.003 


6.45 


91 


-0.21 


0.33 


5.73 


-0.13 


0.20 


5.29 


92 


0.008 


-0.02 


5.92 


-0.01 


0.003 


5.72 


93 


0.005 


-0.02 


5.87 


0.014 


-0.023 


7.83 


fl 


-0.06 


0.023 


5.22 


-0.029 


-0.018 


5.12 




-0.16 


0.29 


6.47 


0.069 


-0.017 


5.69 


fl 


-0.18 


0.25 


8.81 


0.084 


-0.023 


5.13 


9l 


-0.15 


0.14 


5.02 


-0.01 


-0.028 


5.11 


T 
92 


-0.20 


0.31 


5.24 


0.026 


0.04 


4.72 


„T 
93 


0.17 


-0.25 


5.76 


0.11 


-0.18 


5.33 



TABLE VII: Parameters appearing in the fit function of the form factors in full theory for Ej — > 





setl 


sct2 




a 


b 


rrifu 


a 


b 


rrifit 


/l 


0.08 


0.097 


1.53 


-0.12 


0.18 


1.52 


h 


-0.009 


-0.056 


1.57 


-0.01 


-0.029 


1.53 


h 


-0.025 


0.012 


1.61 


0.008 


-0.047 


1.56 


91 


-0.015 


0.31 


1.59 


-0.038 


0.21 


1.60 


92 


-0.008 


-0.12 


1.55 


0.002 


-0.14 


1.61 


93 


-0.026 


-0.13 


1.53 


-0.024 


-0.14 


1.52 


fT 


-0.23 


0.19 


1.52 


0.09 


-0.097 


1.58 


fl 


0.066 


0.067 


1.63 


0.12 


0.13 


1.55 


fl 


0.15 


0.006 


1.56 


0.21 


0.032 


1.65 


T 

9i 


-0.45 


0.29 


1.59 


-0.17 


0.09 


1.62 


T 
92 


0.009 


0.08 


1.57 


-0.026 


0.14 


1.59 


T 
93 


-0.09 


-0.11 


1.54 


-0.07 


-0.16 


1.56 



TABLE VIII: Parameters appearing in the fit function of the form factors in full theory for Ec 
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setl 


set2 




a 


b 


mfu 


a 


b 


nifit 


/l 


-0.22 


0.4 


4.96 


0.037 


0.06 


5.13 


h 


0.009 


-0.024 


5.13 


0.021 


-0.028 


5.28 


h 


0.009 


-0.02 


5.72 


-0.003 


-0.002 


5.67 


91 


-0.029 


0.19 


5.32 


-0.18 


0.28 


5.57 


92 


0.008 


-0.02 


5.41 


-0.01 


0.003 


5.35 


93 


0.005 


-0.018 


4.87 


0.013 


-0.021 


5.06 


fl 


-0.065 


0.025 


5.16 


-0.03 


0.019 


5.25 




-0.24 


0.43 


5.04 


0.104 


-0.026 


5.13 


fl 


-0.19 


0.27 


5.13 


0.091 


-0.025 


5.54 


T 

9i 


-0.14 


0.13 


5.11 


-0.009 


-0.011 


5.14 


T 
92 


-0.03 


0.17 


5.58 


0.04 


0.06 


5.47 


T 
93 


0.21 


-0.27 


5.16 


0.1 


-0.17 


4.96 



TABLE IX: Parameters appearing in the fit function of the form factors at HQET hmit for E;, —¥ nl'^t . 





setl 


sct2 




a 


b 




a 


b 


mfit 


h 


0.02 


0.13 


1.64 


-0.17 


0.25 


1.55 


h 


-0.011 


-0.077 


1.76 


-0.014 


-0.04 


1.51 


h 


-0.034 


0.017 


1.73 


0.011 


-0.065 


1.62 


91 


-0.021 


0.43 


1.68 


-0.053 


0.29 


1.65 


92 


-0.008 


-0.12 


1.57 


-0.002 


-0.14 


1.63 


93 


-0.26 


0.10 


1.62 


-0.022 


-0.13 


1.57 


fT 


-0.097 


0.05 


1.58 


0.098 


-0.1 


1.65 


fl 


0.099 


0.1 


1.69 


0.18 


0.196 


1.59 


fl 


0.14 


0.009 


1.51 


0.32 


-0.048 


1.71 


T 

9i 


-0.4 


0.26 


1.66 


-0.15 


0.08 


1.57 


T 
92 


0.014 


0.12 


1.63 


-0.039 


0.21 


1.63 


T 
93 


-0.08 


-0.1 


1.58 


-0.064 


-0.15 


1.60 



TABLE X: Parameters appearing in the fit function of the form factors at HQET limit for Ec pt^t 
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Full Theory 


HQET 




setl 


set2 


setl 


set2 


/i(0) 


0.14 ±0.04 


0.07 ± 0.02 


0.19 ±0.05 


0.10 ±0.03 


/2(0) 


-0.012 ± 0.003 


-0.006 ± 0.002 


-0.014 ± 0.004 


-0.005 ± 0.002 


/3(0) 


-0.013 ± 0.003 


-0.006 ± 0.002 


-0.011 ± 0.002 


-0.005 ± 0.002 


Si(0) 


0.12 ±0.03 


0.07 ± 0.02 


0.17 ±0.04 


0.10 ±0.03 


52(0) 


-0.012 ± 0.003 


-0.007 ±0.002 


-0.012 ± 0.004 


-0.007 ±0.002 


53(0) 


-0.014 ±0.004 


-0.009 ± 0.003 


-0.013 ±0.004 


-0.008 ±0.003 


/r(o) 


-0.03 ±0.01 


-0.04 ±0.01 


-0.03 ±0.01 


-0.010 ±0.003 




0.13 ±0.04 


0.052 ±0.020 


0.19 ±0.05 


0.079 ±0.003 




0.07 ± 0.02 


0.061 ± 0.020 


0.08 ± 0.03 


0.066 ± 0.021 


alio) 


-0.012 ± 0.003 


-0.03 ±0.01 


-0.012 ± 0.004 


-0.020 ± 0.006 


glio) 


0.11 ±0.03 


0.066 ± 0.021 


0.16 ±0.04 


0.10 ±0.03 


gl{0) 


-0.07 ± 0.02 


-0.073 ± 0.025 


-0.07 ± 0.02 


-0.066 ± 0.021 



TABLE XI: The values of the form factors at = for ne+i'. 





Pull Theory 


HQET 




setl 


set2 


setl 


set2 


/i(0) 


0.19 ±0.05 


0.05 ± 0.02 


0.16 ±0.05 


0.069 ±0.022 


/2(0) 


-0.066 ±0.021 


-0.04 ±0.01 


-0.078 ±0.023 


-0.047 ±0.014 


/3(0) 


-0.013 ± 0.003 


-0.039 ±0.012 


-0.010 ± 0.003 


-0.034 ±0.011 


fli(O) 


0.30 ± 0.09 


0.17 ±0.06 


0.40 ± 0.12 


0.24 ± 0.06 


52(0) 


-0.12 ± 0.03 


-0.14 ±0.04 


-0.12 ± 0.04 


-0.14 ±0.04 


53(0) 


-0.16 ± 0.05 


-0.16 ±0.05 


-0.15 ± 0.05 


-0.15 ±0.04 


/f(0) 


-0.039 ± 0.012 


-0.007 ± 0.002 


-0.042 ± 0.013 


-0.0020 ± 0.0007 


fno) 


0.14 ±0.04 


0.25 ±0.07 


0.21 ±0.06 


0.38 ±0.12 




0.15 ±0.05 


0.24 ±0.07 


0.16 ±0.04 


0.26 ± 0.08 




-0.16 ±0.05 


-0.08 ± 0.03 


-0.14 ±0.05 


-0.07 ±0.02 


glio) 


0.09 ± 0.03 


0.10 ±0.03 


0.14 ±0.05 


0.15 ±0.05 


gim 


-0.20 ± 0.07 


-0.23 ± 0.06 


-0.18 ± 0.05 


-0.21 ±0.08 



TABLE XII: The values of the form factors at = for Sc pe+£~. 
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and 

ri is) 



-2r 



+6mft 



-2mQp ( 1 + 2^ ) ReCf C; +5 



1 + 2^ 



eff 



^io\ 



m. 



4t^ I 1 



-f +3(l + r)i 



r'^ff 



10 



\C- 



101 



T,{s) = 6V^il-t)U[l + 2^\mlp\C, 



-ps 



-12 ( 1 + 2-i ) rriQ [t - r) (1 + sp') ReCf C?, 



9 

9^ 



l-2^)|Cio| 



r3(s) = 12 1 + 2^ rngV^spReCfC; 



2t^ 1 



2^j+4.(^l-^)-3(l + .)t 



sp 



eff 



6mf (2r - (1 + r) i) 



/ 2mQ 



(l^9"l 



Here, G_f 
rriQ 



1.17 X 10 ^ GeV ^ is the Fermi coupling constant, / 



/i+gi „ 

2 ' " 



/i+gi ' 



/i-gi 



(29) 



(30) 



\Cv 



(31) 



2mf 



"["^Sq +'^Ar ~ 9^] ^-nd is the lepton mass. For the Wilson coefficients, 

we use C7 = -0.313, Cg = 4.344, do = -4.669 [32!|. Here we should mention that the Wilson coefficient Cg^^ 
receives long distance contributions from J/ip family, in addition to short distance contributions. In the present 
work, we do not take into account the long distance effects. The elements of the CKM matrix Vtb = (0-77 -0 24)^ 
Vtd = (8.1 ± 0.6) X 10-^ Vfcc = (41.2 ± 1.1) X 10^^ and Vbu = (3.93 ± 0.36 ) 10' ^ have also been used [13]. 

Using the formula for the decay rate the final results as shown in Table IXlIIl are obtained. From this table, we see 





E(, — > ne^e 


Eb — > np,^ fi 


Ef, — > riT^T 


Ec — > pe^ e 


Ec — f pp^p 


Full (setl) 


(4.26 ± 1.27) X 10"^" 


(2.08 ± 0.70) X 10"^° 


(1.0 ±0.3) X 10"^^ 


(5.59 ±1.78) X 10"^^ 


(9.7 ±2.7) X 10"^*' 


Full (set2) 


(5.4 ± 1.6) X 10"^^ 


(2.64 ±0.79) X 10"^^ 


(4.01 ± 1.25) X 10"^^ 


(1.35 ±0.35) X 10"^^ 


(2.36 ± 0.80) X 10"^® 


HQET(setl) 


(8.20 ± 3.04) X 10"^° 


(4.25 ± 2.07) X 10"^° 


(6.26 ±2.46) X 10*^^ 


(7.99 ± 3.07) X 10"^^ 


(1.50 ±0.58) X 10"^^ 


HQET(set2) 


(1.10 ± 0.33) X 10"^" 


(5.67 ± 1.73) X 10"^^ 


(1.16 ±0.46) X 10"^^ 


(2.50 ±0.81) X 10"^^ 


(4.30 ± 1.36) X 10"^" 



TABLE XIII: Values of the r(E(,,c — > ) in GeV for different leptons and two sets of input parameters. 

that: a) The value of the decay rate decreases by increasing in the lepton mass. This is reasonable since the phase 
space in for example r case is smaller than that of the electron and p cases, b) The order of magnitude on decay 
rate of bottom case shows the possibility of the experimental studies on the Sb — > n l^l^ transition, specially the 
p case, at large hadron collider (LHC) in the near future. The lifetime of the Yij, is not exactly known yet but if 
we consider its lifetime approximately the same order of the b-baryon admixture (Ab, S;,, St, ^^b) lifetime, which is 
T = (1.319[j'Q3g^^) X 10""'^^ s [S^, the branching fraction is obtained in 10~^ order. Any measurements in this respect 
and comparison of the results with the predictions of the present work can give essential information about the nature 
of the Sq baryon, nucleon distribution amplitudes and search for the new physics beyond the standard model. 
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Appendix 

In this Appendix, the exphcit expressions for the form factors /i and /2 for b case as well as the nucleon DA's are 
3;iven: 



HQ') 



1 



V2Xi 



1 fl—X2 1 

to "^0 



2\/2 



m6Ul + 3/3)7^19 (xi) -2(-l + ^)Hi7(xi) 



-(3 + P)H5{xi) I - mjva;2 j^i^ -11-13,192, -5,7(2;,) + ^Wii,i3,-i72,-i98,32,5,-7(a;i)| 



+ 



/ dx2 / dxi / dhe-^^'-'Q^y''- 

Jtn Jo Jto 



m%mb 



{-1 + l3)X2n22{Xi) 



m%X2 



M%tl2%/2 
-m^m6(-l + /3)(2 + 3x2) H22{xi) 



{-1+ l3)H-w,iQ{xi) + 2PH2A{xi) + Am%X2+mb{Q'^ + s{ti,Q'')}{-'^ + P)x2 



I'M 



A^4^l2^/2 



TOArX2 



Q2 + s(ti, 16,-224 + (-1 + I3)nw{xi 



-/J'H 16,242 (a^i) +TO6(-l + /3) g2(i + 3a.2) + s(ti,Q2)(l + X2) ■H22(.^;,:) + mlrTOfc ^2(1 + 3/3)7^i6(xi) 



+2(-l + l3)H24{Xi) + (3 + P)Hio{Xi) - (-1 + ^)(3 + X2)U22{Xi) 



m 



N 



{-1 + + X2)Hl0,-16{Xi) 



2{/?(l + X2)n24{x^) + (2 + 4X2)H22(X,)} 



'TV 



+mj^mb 



MpV2 I 

(-1 + /3)W22,-242 (a^O - (1 + 3/3):Wi6(xO - (3 + l3)'Hw{xi) 



SmbQ\-l + |3)'H22ix^) 



{-8 + 2ti-2x2)H22{xi) 



+ (-1 + P)nio,-16{Xi) - 2pn24{Xi) 



+ ruN 



Q2(-1 + /3)(-l + ii - X2)^io,-i6(a;i 



+Q^(-Gti + 6x2 + 4 + 2^)-H22(a;i) + 2Q^P{1 -ti+ X2)H2i{xi) 



Mltl2V2 



U 



6,-183,20 



(Xi) 



+(-1 + p)Hi2{xi) - He,-i3Mxi) 



niN 



M%ti4:V2 



[Q^ + s{ti,Q^)] (3 + 25p)n2o{xi) + 2(-l + p)U-e,i2{xi) 



(5 + /3)Hi8(xi) - m]^ 2(-l + P)n-(i.i2{xi) - (11 + 3/3)7^i8(x,) + (5 + 67/3)%o(x,) 



+2X2 



(-1 + P)U-W,i&{Xi) + m2i{Xi) 



+^'^62,-8,-9,126,143,153,204,-214 (a^i 



M|4\/2 



■^66, -83, -93, 122, 14,15, -204,214 (a^i) + 4X2(-1 + l3)T-L22{Xi) 



^-62,+122,189,-203(a;i) + ^'^62, -122, 183,-2047(2=1 



+4mAr(-l + l3)'H22{xi) + s(ii, Q^) 'Hi83,-2o(a;i) + /3'Hi8,-202i (a^i) 

/3'H44,8,-9,-102, 14,-15,162, -18,2041, -2l4,23i6,244(a;i) + 'H-24, -8,9,102, -14,15,-162, -183,20,-234 (a^i) 



+8{h - X2)H22{Xi) 

rriN 



ti4V2 



2(-l + /3)H-6Mxi) + (1 + 5l3)H2o{xi) - (3 + /3)Wi8(a; 



4 



+ 



4\/2 



(1 + 21p)n2o{xi) - (3 + misixi) 



+ 



/•I /-l- X2 

/ (ix2 / dxie-'"'/^^ 

An Jo 



4 +2 



iQ^ + m%tlfV2 



{to - X2) 



m%mbi-l + /3)(-2 + to){-l + h) + 2m^to{2 + (-4 + to)to} - 2mNtl{Q'^i-2 + 3to) + (-2 + io)s(so, Q^)} 
-mb{-l + /3)to{Q^(-l + 3to) + (-1 + to)s(so, Q^)}] •H22(a;i) + mNto(^{m%{-l + /3)(-l + io) 
- mjvm6io(3 + ^) + (-1 + /3)to[Q^(-l + to) - s(so, Q^)]}^io(a;i) - [m%{-l + /3)(-l + to) + mivm6to(l + 3^) 
+(-1 + p)to{Q'^{-l + to) - s{so, (3^)}1 ^i6(a;i) + 2 [ - niNmbtoi-l + P) + mj^Pil - to) + /3to(Q^(l - ^o) 



m 



JV 



(to - X2) <^ m%mb{-l + /3){-2 + to){-l + to) 



+sM')) n24ix.), , , (^,^^^^,^,272^ 
+2m^to{2 + (-4 + to)to} - 2mNtl{Q^{-2 + 3to) + (-2 + to)s(so, Q^)} - mbto{-l + ^){g2(-l + 3to) 
+(-1 + to)s(so, Q')}] H22(a:0 + mjvto [{m^(-l + /?)(-! + to) - m.Nm.bto{3 + + (-1 + ^)to[Q'(-l + to) 
-s{so, Q'')]}'Hio{xi) - {m%{-l + /3)(-l + to) + mjvm6to(l + 3^) + (-1 + /3)to((5^(-l + to) - s{so, Q''))}Uw{xi) 
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2mN{tQ - X2) m^mbi-l + /3)(-2 + to)(-l + io) + 2mirto(2 + (-4 + io)io) 



(g2+m2,t2)4V2M|to 

+m6to(-l + - 3to) + 2M|io + (1 - to)s(so, g')} + 2mNtl{Q^{2 - 3to)} + 2M|to 

+(2 - to)s(so, Q^) T-L22{xi) + to TOAT^s'Hea, -122, -186,202(2^8) + /3'H-62,i22,-i82,2026(a;i) + 2to;^(1 - to)'Hio,i6(a;i) 

+?Ti^M^to'H-62422, 1811,-205 (2^1) + 'TlArmbM^fo'H-66,83,93,-122,-14,-15,204,-2l4(a;i) + M^Q^to^62,-122,-185,203 (a^i) 
+ mjv;0to'H-io2, 162,244 (x,:)+m^M|;3toH62, -122.183,-2067(2^0 + '7iW'7t(,A^_B/3to'H-62,8,9,-126,-143,-153,-204,2l4(a;i) 
+M|(5^/3ioH-62, 122, -18,2025 (a-'i) + 2TO^io'H-io,i6(a;-j) + m^TO6to'H_io6,-i62,244(a;i) 

+m|fMfto'H_24, -8,9,102, -14,15, -162,-183,20,-234(2:1) + 2m|rQ^to'Hio,-i6(a;i) + Mg(3^fo^_62, 122, 189,-203(2:1) 
+m^/3fo'^i02,-i62,-244(2;i) — ?7iir "^6^*0^102, 166,244(2:1) + ?7iArAf£/3fo'H44,8,-9,-i02, 14,-15, 162,-18,2041, -214,2316,244(2:1) 
+m^Q^/3<o'^-i02, 162,244(2:1) + M|(3^^to'^62, -122, 183,-2047 (2:1) - 2(3^tQ'Hio,-i6(2;i) + Q^/3to'Hio2,-i62,-244(2;i) 
+M|ios(so,Q')H62, -122,-185,203 (a-t) + M%l3tQs{so, <3^)'H_62, 122, -18,2026 (2:0 + 2TO^tos(so, Q^)'Hio,-i6(2:i) 
+M|tos(so, Q^)'Hi83,-2o(2;i) + TO^/3tps(so, Q^)'H-i02,i62,244(2:j) + Mgl3tls{so, Q^)'Hi8,-202i (2:*) 
-2m^a;2 [{(-I + /3)(-l + io) - rnjvm6to(3 + /3) + (-1 + /3)<o[-A^| - Q^(l - to) + s(so, 0^)]}-Hio(2;i) 

- m%{-l + /3)(-l + to) + m;vm6to(l + 3/3) + (-1 + p)to{-Ml - Q^{1 - to) - s{so, Q^)}]nie{xi) 



+2 



- mNmbto{-l + m%l3{l - to) + l3to{Ml + Q^{1 - to) + s{so, Q^)} U2i{xi) 



(32) 



HQ^) 



1 



/ dX2 
J to Jo 



I—X2 



\/2Ae, 

dx2 / dxi / dtie-^(*i'« )/^^ 

to Jo J to 



dxie 



S(X2,Q')/Mi 



2^/2x2 



■^11,-172,5(2:,) - /J'Hii, -1712,5(2:0 





+2(-l + p)H2i{xi) + (3 + mw{xi) 

^,2 



m 



N 



-M + P)X2n22{Xi) + 



m 



N 



m.NmbX2 



+ {b + P)x2 'H22{X,) 
+2 



'■N 



M|^l2^/2 
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2 

mNrnj, 



M%tl%/2^ ' ' ' M%t\2sf2 

- mNmbX2{-l + /?) - {Q^ + s(ti, Q2)}(3 + ^)a;2 + m^(3 + /3 

(1 + 3p)nie{xi) + 2(-l + P)n24xi) + (3 + (3)nw{xi) 



(1 + 3/3)^16(2:0 



mjvmfe(l - /3) - s(ti, Q^)(3 + /? + xa) + m^(5 + /? + xa) - + ^ + (4 + p)x2) H22{x^) 

vv,2 



+ 



"N 



M^V2 

2{2 + j3)x2'H22{x 



- (4 + /3)Q2 - s(ti, 0^) W22(x0 \ + 



M|tf2\/2 



- mb 



'^Wbs, 12, -18, -20(2:1) + /3'H6,123,18,2o(2:i) 



1 Tfi^ r 

J M^t 4v/2 1^-24-8,9,15,-182,-202,228,-234(2:1) + (-1 + l3)'Hu{xi) 



—/3'H-44, -8,9,15,182, -2022,214, -224,-2316,(2:0 ~ 8'H22(2:i) \ H — — "^22(2:, 



^1 I-1-X2 2 

+ / dx2 / dxie-'"'/^^ 

to Jo 



iQ^+m%tl)^V2 



} 

{to - X2) 



m 



V2. 



niNmbto 



(1 + 3^)Wi6(2;i) + 2(-l + l3)H24{xi) 



+ (3 + P)Hw{x^) + 2{mNmb{-l + /3)to + m%{3 + /3 - (5 + p)to + t^) + to{Q'^{3 + /3) - (4 + 0)to) 



+{S + l3-to)s{so,Q^) H22{xi) 



m 



N 



{to - X2) { mNmbto (1 + 3|3)nle{x^) + 2(-l + p)H24{xi) 



+{3 + p)nio{xi 



+<0(-3 -P + to)s{so, Q') U22{Xi) ) + 



{Q-'+mltl)-'2V2' 

- mNmb{-l + P)to + QHo{-3 - /3 + (4 + p)to} + m%{-3 + /3(-l + to) - (-5 + to)to} 

TUN 



(Q2 + m%tl)As/2M%to 



AmN{to - X2) niNmb{-l + P)to + mjv{3 + P 
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-(5 + p)to + 4} + io{M|(2 + /3 + 2io) + Q2{3 + /3 - (4 + /3)to} + (3 + ^ - to)s{so, Q^)}n22{xi) 



-to 



mbMl{ne^^i22,-2Q2{xi) + ^^62,126,202 (a^i)} + 'm%mbto{niOe,i62 -24Axi) + ^^102,166,244 (a;*)} 



,14,-15,202,234(2^1) + /3'^-44, -8,9,-14,15,-2022,214, -23i6 

+miv(l + ^)to}^i8(a;i) - 2m|fm6a;2{(3 + ^)Uw{xi) + (1 + 3/3)Wi6(a;i) + 2(-l + p)U2i{xi)} 



(33) 



where 



H{xi) = H{xi,X2,l — xi - X2) 

-.2X /I .2 , (1-y) 



s(y,(3') 



2/ 



The to = io(sO) Q'^) is the solution of the equation s(to, Q^) = sq, i. 



c., 



(34) 



to(ao,g') = ^ 



+ V4m^(TOg + Q2) + {jnl - - sq^ - sp 
2m% 



(35) 



Here, ,so is continuum threshold, Mg is the Borel mass parameter. In calculations, the following short hand notations 
for the functions %±i^^±j^,... = ±aHi ± bHj... are used, and Hi functions are written in terms of the DA's in the 
following way: 



^1 = ^1 


'H2 = <S'i,-2 


^3=^1 


= Pl,-2 


= Vi 


He = Fi,-2,-3 


n7 = V3 


■Hs = -21^1,-5 + 1^3,4 


% = 14,-3 


Hio = — Vi_-2,-3,-4,-5,6 


^11 = Ai 


H12 = — ^1,-2,3 


= A3 


Hu = -2^i,_5 - ^3,4 


Hi5 = A3,-4 


"Hie = ^1, -2, 3, 4, -5, 6 


Hi7 = Ti 


His = 7i,2 — 2T3 


7^19 = Tr 


7^20 = ^1,-2 — 2T7 


U21 = -Ti,_5 + 2T8 


■H22 = T2, -3, -4, 5, 7,8 


'H23 = ^7,-8 


'^24 = — 7l,-2,-5,6 + 2T7_ 



(36) 



where for each DA's, we also have used X±i^±j^,,, = ±Xi ± Xj.... 
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The explicit expressions for the nucleon DA's is given as: 

V2{xi,n) = 24xiX2[0^(M)+'/>,t(/i)(l-5,T3)], 
V3{xi,n) = 12x3{V'4(^)(l-x3)+-0J(/i)[a;i+X2- 0:3(1 -xa)] 
+■04 (a<)(1 - X3- 10a;ia;2)}, 

V4{Xi,n) = 3{ll)^{n){l-X3)+ll)^{lJ.)[2xiX2-X3{l-Xs)] 
+7P+{f,)[l-X3-2{xl+xl)]}, 

V5 (x,: , fj.) = 6x3 [(i^ (/i) + (/'-) (1 - 2.X3)] , 
Ve{x^,fi) = 2[0O(^) + 0+(^)(l-3x3)], 

Ai(a;j,/Lt) = 120X1X2X303 (/x)(x2 — Xi), 

A2{Xi,n) = 24xiX204 (/i)(X2 - Xi), 

A3{Xi,tx) = 12x3(x2-xi){(V'2(M)+V|(M))+V'r(M)(l-2x3)}, 

^4(Xi,Ai) = 3{x2-xi){-tpl{n) + il)^{fj,)x3+il)^{iJ,){l-2x3)}, 

A^ixi,)!) = 6X3(X2 - Xi)0^(^) 
A6(xi,/^) = 2(X2 - Xi)0g (/x), 

ri(xi,M) = 12OxiX2X3[0g(M) + ^(0^ -0^)(m)(1-3x3)], 
T2{x,,^l) = 24xiX2K2(M)+Ci"(M)(l-5x3)], 

r3(x„ m) = 6x3{(e2 + </'4 + V'4)(m)(1 - + (S4 + -^4 - ^4')(M)[a;? X3(l - X3)] 

+ + <Pt+ - 3^3 - 10X1X2)}, 

r4(x„A*) = ^{(C? + 0^ + V'5°)(m)(1 - ^3) + (^5^ + 4>5 - V'5")(m)[2xiX2 - X3(l - X3)] 

+(4+ + 'Pt + ^Mi^ - ^3 - 2{xl + xi))}, 

T5(Xi,M) = 6X3[e^(/x)+^5+(M)(l-2x3)], 

n{x,,fi) = 2[0°(/i) + i(06--0+)(Ai)(l-3x3)], 

T7(Xi, ^) = 6x3{(-^2 + 04 + V'4)(m)(1 - 2:3) + i-ii + n - ^i ){f^)[xl +xl- X3(l - X3)] 
+ {-Ct + <l>t+ ^4")(m)(1 -X3- 10X1X2)}, 

T8ix^,fi) = ^{(-e? + 0^ + VA°)(A*)(1 - X3) + (-4" + - i^^){l^)[2xiX2 " X3(l " X3)] 

H-U +<Pt+ V'5+)(m)(1 - ^3 - 2{xl + xl))}, 
Slixi, n) = 6X3(X2 - Xl) [{^l +ct>l + i)l+^t + 4^ + i't){lJ) + (C +4i- V-DIm)!! - 2x3)] 

52(x„ A*) = ^(X2 - Xl) [- (^Ag + 00 + it) (m) + (4" + <P5 - V's) (M)a:3 

+ (e^ + 0^+V5°) (m)(1-2x3)] 

Pl(x„ m) = 6x3(x2 - Xl) [(^2 -41-4^1+ -<t't- V'DCm) + (C - <^4 + V'4")(m)(1 - 2X3)] 
P2{XU m) = ^(X2 - Xl) [{^l + Vs" - Cs) (m) + fe" - '/'S" + V-s) (M)a:3 

+ fe+-0+-^°)(M)(l-2x3)] . (37) 
The following functions are encountered to the above amplitudes and they can be defined in terms of the eight 



independent parameters, namely Jn, Ai, A2, Vf, A\, /J and 







^6 = /JV 


<t>l 


= 


</'^ = ^(Ai + /iv) 


a 


= 


4 - gA2 




= 


V'5° = ^(/iv-Ai) 




= 






= 




<i>4 


= 


^(Ai(l-2/{'-4/f) + /;v(2A^-l)) , 


<i>t 


= 


J(Ai(3-10/|')-/;vW-3)) , 




= 


-J (Ai (2 - 7/1^ + /r) + /^(A^ + - 


i^t 


= 


-i (Ai(-2 + 5/1'' + 5/1") + /jv(2 + 




= 


Aa2(4-15/|), 




= 


^A2(4-15/2''), 




= 


^(Ai(/f-/r) + /iv(2^5'-l)) , 




= 


-^(Ai(4/f-l) + M3 + 4y/')) , 


V'5- 


= 


^(Ai(A'-/r) + /^(2-A^-3yi'^)) , 




= 


(Ai(-1 + 2f^ + 2/n + /iv(5 + 2^« 




= 




4+ 




|a2(2-9/,''), 






i(Ai(l-4/f-2A«) + /;v(l+4A^)) 






-i(Ai(l-2A^) + M4Fi'^-l)) 



